function [result] = read3(input,dim1,dim2,dim3)
    result(dim1,dim2,dim3)=0.0;
    a=dim3*dim2;
    b=dim3;
    for h=1:dim1
        for i=1:dim2
            for j=1:dim3
                result(h,i,j)=input(a*(h-1)+b*(i-1)+(j));
            end
        end
    end
end

